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Abstract 


This  report  documents  our  procedures  and  the  accuracy  of  objective  maps  of  bottom 
pressure  in  the  Synop  Central  Array  region  between  June  1988  and  August  1990.  The 
pressure  maps  are  made  by  using  multivariate  objective  analysis  procedures  to  combine 
data  from  bottom  pressure  measurements  by  Inverted  Echo  Sounders  (PIES)  and  bottom 
current  meter  measurements  in  the  Central  Array  between  June  1988  and  August  1990. 

A  key  new  feature  is  that  all  daily  maps  can  be  adjusted  to  a  fixed  reference  level  by 
combining  measurements  of  bottom  pressure  and  near-bottom  currents,  (i\,  tt,v).  While 
objective  stream  functions  mapped  from  (o,v)  alone  are  spatially  consistent,  the  reference 
levels  vary  temporally  from  map  to  map.  On  the  other  hand,  each  bottom  pressure  record 
Ph  has  a  temporally  consistent  reference  but  the  reference  can  vary  from  site  to  site.  The 
combination  allows  us  to  adjust  all  to  a  common,  consistent  reference  level  in  multivariate 
objective  maps. 

The  pressure  maps  and  error  maps  are  displayed  daily  for  26  months  with  overlaid 
measured  current  vectors.  Mapping  procedures  are  documented  along  with  error  analyses 
and  comparisons  with  measured  fields.  The  results  produce  bottom  pressure  and  current 
fields  with  typical  error  of  only  2  mbar  and  2  cm/sec,  compared  to  typical  signal  standard 
deviation  of  6  mbar  and  9  cm/s. 

The  report  is  one  of  a  series  of  data  reports  that  document  aspects  of  our  participation 
in  the  SYNoptic  Ocean  Prediction  (SYNOP)  experiment. 
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1  Introduction 

This  report  documents  the  procedures  &nd  the  accuracy  of  the  objective  bottom  pressure 
maps  in  the  Synop  Central  Array  region  between  June  1988  and  August  1990.  It  is  one  of 
a  series  of  data  reports  that  document  aspects  of  our  participation  in  the  SYNoptic  Ocean 
Prediction  (SYNOP)  experiment  [7,  2,  3,  8,  9,  6].  The  study  region  consists  of  a  sub- region 
of  the  Synop  Central  Array  and  is  shown  in  Figure  1  and  Table  1. 

The  bottom  pressure  maps  are  made  by  combining  data  from  bottom  pressure  measure¬ 
ments  [2,  3]  by  Inverted  Echo  Sounders  (PIES)  and  bottom  current  meter  measurements. 
The  instrument  positions  are  indicated  in  Figure  1  and  detailed  in  Table  2  and  Table  3.  For 
each  instrument  site,  the  period  when  PIES  and  Current  Meter  measurements  were  taken  is 
shown  by  the  time  lines  in  Figure  2  and  by  the  start  and  end  dates  in  Table  2  and  Table  3. 

The  rest  of  the  report  is  organized  as  follows.  Section  2  reviews  the  aspects  of  multi¬ 
variate  Objective  Analysis  procedures  necessary  for  mapping  the  pressure  fields  from  the 
pressure  and  current  measurements.  Section  3  describes  the  data  processing  steps.  Sec¬ 
tion  4  evaluates  the  precision  of  the  resulting  pressure  maps.  The  pressure  maps  are  then 
displayed  in  section  5  along  with  the  corresponding  error  maps. 

Table  1:  Corner  positions  of  the  rectangular  mapping  region 


Corner 

X  (km) 

y  (km) 

Latitude  (N) 

Longitude  (W) 

Upper  Left 

-160 

38"  40.27 

70"  09.10 

Lower  Left 

-160 

36"  14.08 

69"  18.69 

Upper  Right 

40 

120 

39"  08.25 

67"  54.73 

Lower  Right 

40 

36"  42.06 

67"  04.32 

2 


Figure  1:  Current  meter  mooring  and  PIES  sites.  The  dashed  curve  indicates  the 
path  of  the  Gulf  Stream  (1975  to  1986)  from  Gilman  and  Comillon  [1990]. 
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Figure  2:  Time  lines  of  the  data  coverage  for  current  meters  and  PIES’s  on  all  sites 
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Table  2:  Duration  and  location  of  the  PIES  records  used  in  mapping 


Site 

Start  date 

End  date 

Latitude  (N) 

Longitude  (W) 

Depth  (m) 

G2 

28-NAY- 1988 

12-SEP-1989 

37“  47.52 

69“  24.33 

4060 

07-JUM-1989 

19-AUG-1990 

37“  47.84 

69“  24.31 

4088 

G3 

28-MAY- 1988 

28-MAY- 1989 

37“  17.40 

69“  14.53 

4358 

30-MAY- 1989 

18-AUG-1990 

37“  16.99 

69“  14.71 

4358 

H2 

28-MAY- 1988 

Ol-JUM-1989 

38“  37.88 

68“  54.20 

3458 

17-AUG-1989 

10- AUG- 1990 

38“  37.78 

68“  54.90 

3458 

H3 

13-JUM-1988 

14-AUG-1989 

38“  10.10 

68“  43.30 

4030 

17-AUG-1989 

14- AUG- 1990 

38“  10.09 

68“  43.65 

4025 

H4 

16-Jinf-1988 

14-AUG-1989 

37“  39.80 

68“  35.40 

4445 

17-AUG-1989 

17-AUG-1990 

37“  39.57 

68“  35.35 

4445 

H5 

09-JUN-1988 

13-AUG-1989 

37“  10.60 

68“  17.00 

4790 

16- AUG- 1989 

17-AUG-1990 

37“  10.23 

68“  17.83 

4800 

H6 

24-MAY- 1988 

03-JUN-1989 

36“  40.45 

68“  15.64 

4883 

OS-JUM-1989 

08-AUG-1990 

36“  39.35 

68“  15.70 

4840 

11 

26-MAY- 1988 

Ol-FEB-1989 

38“  47.54 

68“  06.38 

3828 

Ol-JUM-1989 

16- AUG- 1990 

38“  47.58 

68“  06.25 

3828 

12 

ll-JUH-1988 

29-MAY- 1989 

38“  20.90 

67“  59.60 

02-JUM-1989 

15- AUG- 1990 

38“  19.68 

67“  58.71 

13 

07-JUM-1989 

10- AUG- 1990 

37“  47.61 

67“  58.85 

4610 

14 

lO-JUM-1988 

15-JUM-1989 

37“  18.50 

67“  39.30 

18-JUI-1989 

16- AUG- 1990 

37“  18.88 

67“  39.58 

15 

25-MAY- 1988 

02-JUN-1989 

36“  49.73 

67“  27.57 

05-JUI-1989 

07-AUG-1990 

36“  50.19 

67“  27.36 
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T^ble  3:  Duration  and  location  of  the  current  meter  records  used  in  mapping 


Site 

Start  date 

End  date 

Latitude  (N) 

Longitude  (W) 

Depth  (m)  1 

G2 

28-MAY" 1988 

03-JUM-1989 

37®  47.30 

69®  24.50 

06-JUN-1989 

24- AUG- 1990 

37®  47.03 

69®  24.66 

G3 

27-MAY- 1988 

ll-QCT-1988 

37®  17.70 

69®  14.50 

3500 

H2 

26-MAY- 1988 

10- AUG- 1990 

38®  38.02 

68®  54.00 

3355 

H3 

03-JUN-1989 

14- AUG- 1990 

38®  09.31 

68®  44.28 

3500 

H4 

lS-JUN-1988 

OS-JUN-1989 

37®  40.10 

68®  35.20 

3500 

H5 

09-JUN-1988 

14-JUM-1989 

37®  11.05 

68®  17.60 

17-JUM-1989 

23- AUG- 1990 

37®  11.50 

68®  15.96 

H6 

23-MAY- 1988 

2S-MAY-1989 

36®  40.40 

68®  15.85 

3500 

11 

26-MAY- 1988 

lO-AUG-1990 

38®  47.83 

68®  06.25 

3500 

12 

lO-JUM-1988 

29-MAY- 1989 

38®  19.40 

67®  59.20 

3500 

02-JUN-1989 

16-AUG-1990 

00 

o 

00 

67®  59.37 

3500 

13 

lO-JUH-1988 

03-MOV-1988 

37®  47.35 

67®  56.20 

3500 

26-AUG-1989 

09- AUG- 1990 

37®  47.81 

67®  57.82 

3500 

14 

09-JUN-1988 

04-SEP-1988 

37®  18.50 

67®  38.30 

3500 

27- AUG- 1989 

lS-AUG-1990 

37®  19.07 

67®  39.58 

3500 

15 

24-MAY- 1988 

07-AUG-1990 

36®  50.10 

67®  27.40 

3500 

M13 

28-AUG-1989 

17-AUG-1990 

37®  29.12 

68®  07.18 

2  Review  of  Objective  Analysis  Procedures 


This  section  reviews  Objective  Analysis  (OA)  procedures.  We  first  review  the  Ganss-Markov 
problem  of  estimating  a  random  variable  from  some  other  random  variables  that  are  cor¬ 
related  to  it  [5].  The  solution  to  this  problem  is  then  applied  to  mapping  the  pressure 
field  from  pressure  and  current  measurements  [1].  Various  computational  aspects  of  the 
procedures  are  discussed  next.  Section  2.3  examines  the  error  model  underlying  the  OA 
procedure  we  used,  followed  by  a  discussion  ot  normalization  of  the  covariance  functions 
in  Section  2.4.  Programs  are  discussed  in  Section  2.5  and  listed  in  Appendix  A.  For  use 
in  later  sections,  some  theoretical  relationships  between  spatial  cross  covariance  functions 
between  pressure  and  currents  are  tabulated  in  section  2.6,  and  their  derivation  is  provided 
in  Appendix  B. 

2.1  OA  basics 
Given  N  random  variables 

^  (1) 
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where  the  superscript  T  denotes  t’:e  matrix  transpose,  \ve  want  to  estimate  another  random 
variable  0  by  a  linear  relation  Between  6  and 

N 

(2) 

«=I 

where  j  -  (hi, 62?  is  selected  in  such  a  way  that  the  resulting  mean  square  error 


P  =  E{(e-e)^} 


(3) 


is  minimum.  Note  we  implicitly  assume  that  E{6}  =  E{4>,)  =  0,  because  otherwise  the 
estimator  in  (2)  can  be  readily  improved  by  adding  a  constant  term  to  it. 

We  now  show  that  the  optimal  estimator  and  the  minimal  mean  square  error  in  the 
above  problem  are  given  by 

e  =  c>i-V  (4) 

P=  (5) 

where  A  and  C  are  respectively  the  covariance  matrix  between  the  (^'s  and  the  cross* 
covariance  matrix  between  9  and  <t/a.  Specifically,  /t  is  an  iV  by  symmetric  matrix 

A  -  E{H>^}  or  Atr  =  E{<t>t4>r}  (6) 


and  C  is  a  row  vector  of  length  N 

C  =  or  Cr  =  E{9<t>r}.  (7) 

To  prove  the  above  result,  we  note  that  the  meam  square  error  P  in  (3)  is  a  function 
of  b,  and  is  minimum  if  all  its  partial  derivatives  with  respect  to  b,  are  zero.  That  is, 
2E{(9  -  9)^}  =  2£{(^  -  =  0,  r  =  1,2,  ...,1V.  To  write  in  matrix  form,  we  have 

E{{9-9)<t>'^)  =  Q  (8) 

Inserting  (2)  into  (8),  we  obtain  =  E{9<p^),  or 

b^A  =  C  (9) 

By  inserting  the  solution  of  (9)  for  b^  into  (2)  we  get  the  optimal  estimator  in  (4).  Next, 
to  obtain  the  least  mean  square  error  in  (5),  we  use  (3),  (2)  and  (8)  to  obtaun 


P  =  £{{9  -  9)9'^}  -  E{{9  -  =  £{(»  -  9)9}  *  £{#^}  -  E{99) 
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Now  using  right  hand  side  of  (4)  to  replace  B  yields 


which  is  the  same  as  (5)  by  the  definition  of  C  in  (7). 

2.2  Extending  basic  formulas 

We  now  generalize  (4)  and  (5)  for  computational  convenience.  Our  main  interest  is  to 
use  the  result  to  estimate  a  physical  parameter  6  at  some  locations,  where  possibly  no 
direct  measurement  is  available,  from  measurements  <i>  of  possibly  different  parameters  at 
instrument  sites. 

First  we  note  that  can  be  measurements  of  different  kinds.  For  example,  we  can 
have  N  measurement  sites  measuring  M  different  parameters  with  Nm  sites  measuring 
parameter  (thus  N  =  53^=1  Nm)-  In  this  report  we  are  interested  in  estimating  pressure, 
B,  at  an  arbitrary  point  from  measurement  of  pressures  and  currents  at  instrument  sites 
(M  =  3).  Thus  an  example  with  Af  =  3  and  N  —  1  can  be 

4>  =  (Pl,P2,P3.«l,nj,Vi,V2)^. 


Secondly,  when  we  need  estimates  of  B  on  several,  say  G,  locations,  we  can  apply  (4) 
and  (5)  to  all  G  interpolation  points  one  by  one 


B,  =  CgA-^<t> 


P,  =  ElSl)  -  C,A-'Cj 

where  subscripts  of  B  denote  interpolation  points,  Cg  denotes  row  vector  C  corresponding 
to  interpolation  point.  Stacking  those  equations  together,  we  can  put  them  in  matrix 
form 


B  =  CA-^(t> 

(10) 

P  =  diagonal[£{^fl^}  -  Gi4“^C^], 

(11) 

with  B  and  P  now  being  column  vectors  of  length  G,  and  C  being  a  G  by  iV  matrix.  For 
example,  for  G  =  4  we  have 


1 

(  E{BM  ...  E{BiM  \ 

B2 

and  C  = 

E{02^}  ...  E{B24>n} 

h 

£{^3^}  ... 

\  9* 

^  £{^4^}  ...  E{94<^s}  > 

(12) 
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Finally,  if  we  want  to  estimate  a  time  series  (say,  of  length  T)  of  9  from  a  corresponding 
series  of  then  under  stationarity  assumption  about  both  9  and  tj>  fields,  covariance  matrices 
A  and  C  do  not  change  with  time  and  we  can  apply  (10)  at  each  time  point.  The  result  can 
again  be  expressed  in  compact  form  as  (10)  with  matrix  sizes  of  9  and  now  being  G  by 
T  and  N  hy  T  respectively.  Equation  (11)  holds  without  any  change  in  this  case  because 
error  estimate  does  not  change  with  time  due  to  stationarity,  as  can  also  be  seen  directly 
from  the  right  hand  side  of  (11).  For  an  example  with  T  =  3,  we  have 

(<9(l),fi(2),^(3))  =  (13) 


where  C  and  each  9(t)  can  be  of  the  forms  in  (12)  and  each  (t>{t)  can  be  a  column  vector  of 
measurements  at  all  instrument  sites. 

To  summarize  the  above  extensions'  we  repeat  here  the  main  formulas  after  the  exten¬ 
sions 


9  =  CA~^<f>  or 

P  =  diagonal[E{d(t)^(0^}  -  CA~^C^]  or 
A  =  E{4>{t)4>{tf)  or 

C  =  £{d(()^0^}  O’- 


S  /  N  \ 

^gt  —  53  ^9^  (  Xrf  )  ’ 

r=l  \#=sl  / 


N 


Pg  —  E{9^{t))  -  ^2  ('gr^rgCgti 


r,$=l 


(14) 

(15) 


A,.  =  E{<l>,{t)Mt)h  (16) 

Cgr  =  £{f>,(t)-^r(0};  (17) 


where  the  dimensions  of  all  matrices  are  as  follows 
9:  GbyT; 

P  :  G  by  1  (It  would  be  G  by  T  if  either  ^  or  0  were  non-stationary); 
A  :  N  by  N; 

C  :  GbyN; 

4>:  N  by  T,  (notation:  =  ^*(0  and  <^(0  >s  the  t*'*  column  in  (f>)\ 

9  :  G  by  T  (notation:  9tt  =  9,{t)  and  9{t)  is  the  t*'*  column  in  9). 


Here  is  an  example  of  (14)  where  N  =  7,M  =  3,  G  =  i,T  =  3: 


/«,(!)  tf,(2)  ^,(3)\ 
tfj(l)  ^2(2)  ^2(3) 

Ml)  M2)  M3) 
\Mi)  M2)  to/ 


E[9xp,]  E[9Mn\  £ltfiP3l 
E[9gpx]  E\9m]  E\971^] 
E[9zpx\  E\9:Mn]  £[«3P3] 
E[9ipi]  Ei\9iPi]  £[to 


£[tfltt2]  ^ 

£[02^3]  £[^t>i]  £[®lU2i 
£[^3«3l  ^'[^1*^2] 

£[tf4U2l  E[94Vi]  £(tfiW2i  > 


’The  exteuioi  can  also  incinde  the  caae  where  loine  d(ld:*),  for  «  >  1,  are  also  naed  in  eeiimatiag  t(t). 
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E\p\P7\ 

E\p\Pi\ 

£[piui] 

£[PiU2] 

£[pit;i] 

E\P\V2]  '' 

-1 

^  Pl(l) 

pi(2) 

Pi(3) 

E\p2Pi] 

E\piP2] 

E\p2Pz] 

£[P2Ui] 

£(p2U2j 

E{p2Vi] 

E\p2V2] 

P2(l) 

P2(2) 

P2(3) 

^[PaPzl 
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2.3  Modeling  measurement  errors 

The  measurement  error  is  assumed  to  be  additive.  Thus  a  measurement  0,  is  the  sum  of 
true  signal  0o<  and  a  noise  term  c,,  assumed  to  be  uncorrelated  to  all  the  signals  in  <t>o, 


<t>t  =  <i>oi  +  f*. 


Also  by  assumption,  from  different  measurement  sites  are  not  correlated;  they  all  have 
zero  mean  but  may  have  different  variances. 

From  those  assumptions  about  measurement  error,  the  covariance  matrices  A  and  C 
can  be  written  as 

A  —  Aq  +  E  (18) 

C  =  Co,  (19) 

where  Aq  and  Co  are  corresponding  covariance  matrices  in  the  absence  of  noise  and  £  is  a 
diagonal  matrix  whose  diagonal  element  is  the  variance  of  e,. 

In  contrast  to  the  measurement  error,  we  will  call  the  least  mean  squaire  error  in  (5) 
mapping  error.  It  is  sometimes  useful  to  check  the  mapping  error  at  a  measurement  site. 
To  this  end,  we  suppose  that  one  of  the  measurement,  say  <t>k,  happens  to  be  the  sum  of  6 
and  a  noise  term  t 

=  ^  +  €, 

where  t  has  zero  mean  and  variance  <r^  and  is  independent  of  B  and  all  other  then 
equation  (5)  reduces  to 

P  =  a\l-<7^A;li), 

where  Aj^jJ  is  the  kth  diagonal  element  of  A~'.  Thus  the  mapping  error  is  smaller  than  the 
measurement  error  at  a  measurement  site  due  to  the  contribution  of  measurements  from 
other  sites. 


2.4  Normalization  of  covariance  matrices 

Computationally,  it  is  convenient  to  normalize  the  covariance  matrices  A  and  C.  We  will 
denote  normalization  by  a  tilde.  Let 


^ 

(20) 

A  =  n;^an;^ 

(21) 

C  = 

(22) 

P  =  Ng-^P 

(23) 

where  N^,  is  a  diagonal  matrix  whose  diagonal  element,  (iV^)«<,  is  the  normalization 

factor  for  ^«(t)  and  Ng  is  a  diagonal  matrix  whose  diagonal  element,  (Ng)gg,  is  the 

normalization  factor  for  dg(t).  Specifically,  we  have 

Mt)  =  6g{t)  =  eg{t)/{Ng)gg 

Ar.  - 

c.  =  ew)M»)  ^ 

Solving  (21)  and  (22)  for  A  and  C  and  then  putting  them  into  (14)  and  (15),  we  get 

0  =  NgCA-^N;'(i> 

(24) 

P  =  diagonal[£{tf(0^(0’’)  -  NgCA-^C'^Ng]. 

(25) 

Those  two  equations  are  often  written  in  more  symmetrical  form  as 

‘0  =  CA-^^ 

(26) 

P  =  diagonal[£{«(t)«(t)^}  -  CA-’C’’] 

(27) 

where  (20)  and  (23)  have  been  used.  In  particular,  if  we  choose  {Ng)gg 

(27)  becomes 

P  =  diagonal[/G)<C  -  CA~*C^] 

=  y/E[0j{t)],  then 

(28) 

Moreover,  if  all  have  same  value  <t^  and  we  choose 


JV^  =  alsxN  Nt  -  (tIoxG’ 


(29) 
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with  I  being  the  identity  matrix,  then  (24)  can  be  written  as 

0  =  Ci'V  (30) 

Note  in  (30)  neither  0  nor  4>  is  normalized  and  in  (28)  P  is  the  mapping  error  as  a  fraction 
of  the  signal  variance  E{0g(t)^}. 

In  practice,  the  normalization  factor  (N^)^  is  often  chosen  to  be  proportional  to  either 
the  total  variance  or  the  signal  variance  E{(f>Q^}.  The  first  choice  appears  to  be 

more  natural  at  first  since  the  total  variance  is  often  more  readily  available  than  the  signal 
variance.  However,  we  have  adopted  the  second  choice  because  it  facilitates  handling  the 
case  where  measurement  errors  differ  at  different  instrument  sites.  Under  this  choice,  the 
measurement  errors  remain  to  be  constrained  in  the  diagonal  elements  of  matrix  A  (refer 
equations  (18)  and  (19)).  Thus  if  different  error  estimates  are  used  for  different  measurement 
sites,  the  second  approach  requires  only  the  diagonal  elements  of  matrix  A  to  be  changed 
while  the  first  approach  necessitates  changes  of  whole  matrix  C  and  all  off>diagonal  elements 
of  matrix  A. 

In  addition  to  the  computational  convenience,  the  normalization  also  simplifies  handling 
the  case  of  non-uniform  variance  in  the  measured  fields.  In  this  case,  one  normalizes  <l>  and  0 
such  that  ^  and  0  fields  in  (20)  have  uniform  variance,  thus  permitting  use  of  homogeneous 
covariance  functions  in  A  and  C  in  (26).  This  is  equivalent  to  using  (14)  with  A  and  C 
calculated  from  some  covariance  functions  whose  values  at  zero  lag  (variances)  are  location 
dependent. 

2.5  Program  prototypes 

Five  programs  are  listed  in  Appendix  A.  Three  of  them  are  for  mapping  pressure  fields 
respectively  from  pressure  measurements,  current  measurements  and  the  combination  of 
pressure  and  current  measurements.  The  other  two  are  for  mapping  current  fields  respec¬ 
tively  from  current  measu’^ements  and  the  combination  of  pressure  and  current  measure¬ 
ments.  The  programs  closely  follow  the  formulation  in  section  2.2.  The  output  error  fidds 
are  expressed  in  terms  of  variance  percentage  of  signal  fields  and  are  calculated  using  (27) 
(after  factoring  out  the  first  term).  The  relationship  between  input  and  output  fields  can  be 
interpretated  according  to  either  (14),  (26)  or  (30)  depending  on  the  nonnalixation  factors 
used. 


A  dose  match  between  the  statements  in  the  programs  and  the  formtiias  in  section  2  2 
makes  the  programs  more  readable,  thus  reduces  the  chance  of  coding  error,  and  facilitates 
their  adaptation  for  new  situations.  We  used  these  programs  as  the  basis  for  more  ^user* 
friendly’  programs  and  to  test  different  correlation  functions.  However,  Uttle  emphasis  has 
been  given  in  these  programs  to  the  convenience  of  input  format,  input  validation  and 
CPU/storage  efficiency.  Also  ‘bogus  data’  is  not  allowed  in  the  input  matrices  to  these 
programs. 

2.6  Assumed  cross-covariance  functions  between  pressure  and  current 

In  order  to  use  the  OA  formulas  (14)-(17),  we  need  to  supply  the  covariance  matrices  A 
and  C.  In  our  application,  the  elements  of  these  matrices  are  calculated  from  the  spatial 
cross  covariance  functions  between  current  and  pressure.  This  subsection  examines  some 
theoretical  relationships  between  these  cross-covariance  functions  for  use  in  section  3.1  to 
obtain  matrices  A  and  C. 

We  assume  that  the  current  is  geostrophic  and  has  little  vertical  shear  between  the 
current  measurement  level  at  3500  m  depth  and  the  pressure  measurement  level  at  the  sea 
floor.  Also  we  assume  that  pressure  field  statistics  are  stationary  in  time,  and  homogeneous 
and  isotropic  in  space,  and  its  spatial  covariance  function  is  a  Gaussian  function.  Based 
on  these,  we  examine  here  the  required  forms  of  the  cross-covariance  functions  between 
current  components  (u,v)  and  pressure  p.  Further  discussions  on  the  above  assumptions 
can  be  found  in  Section  3.1. 

Specifically,  assuming  geostrophy,  the  current  and  pressure  fields  are  related  by 

1  1  ^P 

“  fopdy'  ^  fopdx' 

where  /o(=  2il  •  sin  a)  is  the  Coriolis  parameter  and  p  is  the  density  of  the  sea  water.  The 
values  of  /o  and  p  we  used  correspond  to  the  values  for  the  Central  Array  around  3500  m 
depth: 

/o  =  2fi  sina  =  2  X  7.292  X  10"*  X  sin 38®  =  8.9788  X  10-®  (s-*); 
p  =  1050  (kg/m^). 

The  current  and  pressure  were  measured,  respectively,  at  3500  m  depth  (above  the  bottom 
frictional  layer)  and  near  ocean  bottom.  Thus  with  (31),  we  have  implicitly  used  the  as¬ 
sumption  that  the  current  has  little  vertical  shear  in  the  region  outside  the  bottom  frictional 
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layer  and  between  the  two  measurement  levels.  With  this  assumption,  the  bottom  pressure 
field  will  be  the  same  as  the  pressure  field  at  the  3500  m  horizon  within  an  overall  constant 
independent  of  location  and  time  (see  also  section  3.3.2),  because  a  frictional  boundary 
layer  is  not  expected  to  alter  the  pressure  from  its  geostrophic  value  outside  the  layer. 

From  geostrophic  relation  (31),  we  see  that  the  velocity  field  is  horizontally  non-divergent. 
Thus  a  stream  function  exists  such  that 


(32) 


From  this,  we  can  derive  all  the  cross-covariance  functions  between  stream  function  and 
current  from  the  covariance  function  of  the  stream  function.  The  derivation  is  given  in  the 
Appendix  B  and  the  results  are  listed  in  Table  4. 


Table  4:  Cross-covariance  functions  of  stream  function  and  current.  7,  Og  and  Oy  are  defined 


respectively  by  7  =  {(xj  -  xj)^  -|-  (yi  -  xj  -  *1  =  7cosa*  and  ya  -  yi  =  7co8a,. 

V'(l2,y2)  «(X2,y2) 

»(*a.Va) 

«(xi,yi) 

u(xi,yi) 

^(7)  -COSOj,^ 

-cos2aj,0  -  (l-COS^Oy); 
-cosa*^  COSOrCOSOy  (0  -  ig] 

cosa.^ 

cos  a*  cos  tty  (^  -  ig) 

1  -cos2a*0 -(1 -co8^o*)ig 

Table  5:  Same  as  the  last  table  except  that  F  is  assumed  to  be  a  Gaussian  function  here. 

ti'(*2,y2) 

«(*2,y2) 

»(*2,ya) 

■ 

f(7)  =  V^exp(-A7*)  7cosay  •2Af’(7)  -7CO8Q,  •  2AF’(7) 

-7Cosay  •  2AF(7)  (1  -  2A7*cos*  Oy)  •  2AF(7)  2A7*cosa,co8ay  •  2AF(7) 

7cosai  •  2AF(7)  2A7*cosoiCosay  •2AF(7)  (1  -  2A7*cos*q*)  •  2AF(7) 

From  (31)  and  (32),  we  see  that  the  stream  function  0  can  be  chosen  in  such  a  way  that 
pressure  p  is  proportional  to  V't  namely, 

P  =  /oPV’-  (33) 

Therefore,  all  the  cross-covariance  functions  between  pressure  and  current  are  the  same  as 
those  between  stream  function  and  current  except  for  a  constant  factor  and  thus  can  be 
readily  obtained  from  Table  4. 

In  particular,  we  will  assume  that  the  pressure  covariance  function  is  a  Gaussian  function 


^(7)  =  E{\pi  -  Eipi)]\pj  -  £(Pi)l)  =  Vpexp(-A7") 


(34) 
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where  7  =  -  *j)^  +  (Vi  -  VjY  **  the  distance  between  site  *  and  site  j,  and  A  = 

(Correlation  Length')'*  ^  constant  to  be  fitted  by  the  data,  and  Vp  is  the  signal  variance. 
By  Equation  (33),  the  covariance  function  of  the  stream  function  0  is  also  a  Gaussian 
function.  As  shown  in  Appendix  B,  the  resulting  cross  covariance  functions  between  current 
and  stream  function  will  be  of  the  forms  listed  in  Table  5. 
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3  Mapping  the  pressure  fields 

This  section  describes  the  data  processing  procedures  involved  in  mapping  pressure  fields. 
The  major  work  of  objective  mapping  consists  in  estimating  the  noise  variance  in  measure¬ 
ments,  modeling  the  covariance  functions  and  estimating  the  mean  fields  based  on  measured 
fields.  We  discuss  these  problems  in  the  following  two  subsections.  Then  the  final  subsec¬ 
tion  describes  the  processing  procedures  that  were  applied  to  the  measured  pressures  and 
currents  and  eventually  led  to  a  time  series  of  pressure  maps. 

3.1  Estimating  covariance  functions 

The  noise  variance  in  both  pressure  and  current  measurements  is  chosen  to  be  ten  percent  of 
the  signal  variance  (Table  6)  by  the  following  considerations.  The  typical  standard  deviation 
of  ail  the  pressure  records  is  about  0.06  decibar.  The  error  in  residual  pressure  is  negligible 
compared  to  the  error  in  the  reference  pressure  (see  section  3.3.2).  The  error  standard 
deviation  of  the  reference  pressure  is  estimated  to  be  0.02  decibar.  Thus  for  pressure, 
Vn  «  (o!§i)^  **  0-10.  The  typical  standard  deviation  of  the  deep  current  records  is  about 
0.1  meter  per  second  while  the  error  standard  deviation  is  estimated  to  be  0.03  meter  per 
second.  Thus  for  current  record,  Vn  as  (  o^) 


Table  6:  Noise  variance  and  correlation  length  used 


Fields 

Noise  Variance/ Signal  Variance  Correlation  length 

Pressure  or  current 
Mean  pressure  or  current 

0.10  90  km 

0.05  150  km 

The  cross-covariance  functions  between  pressures  (or  stream  functions)  and  currents 
are  assumed  to  be  of  the  forms  listed  in  Table  5.  In  doing  so,  we  have  made  foUowing 
assumptions  as  discussed  in  section  2.6.  These  assumptions  are  only  approximations  to 
reality  and  they  are  made  to  facilitate  the  data  processing  and/or  make  full  use  of  existing 
measurements  in  improving  the  mapped  pressure  fields. 

•  geostrophy.  The  current  measurements  are  typically  on  the  order  of  0.2  m/sec  with 
a  length  scale  of  about  100  km.  Combining  these  values  with  the  value  of  Coriolis 
parameter  given  below  Equation  (31)  gives  a  Rossby  number  of  about  0.02.  Thus,  we 
expect  that  geostrophy  is  a  good  approximation. 

•  no  vertical  shear  in  currents  between  3500m  depth  where  currents  were  measured  and 
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top  of  the  bottom  friction  layer;  no  change  in  geostrophic  pressure  through  the  bottom 
friction  layer  to  the  ocean  bottom.  Hence  the  PIES  measurement  of  pressure  at  the 
bottom  is  the  same  as  at  the  3500  m  horizon  except  for  a  dynamically  unimportant 
constant  hydrostatic  reference  offset. 

•  homogeneity.  From  our  measurements,  the  fields  on  the  northern  part  of  the  array 
typically  have  less  variance  but  more  high  frequency  components  than  their  southern 
counterparts.  Although  nonhomogeneity  of  other  spectral  characteristics  often  has  no 
simple  solution,  variation  in  variance  in  the  measured  fields  can  usually  be  remedied 
by  normalizing  the  measurement  fields  by  their  standard  deviations,  then  mapping 
the  normalized  fields  using  a  homogeneous  covariance  functions  and  finally  scaling 
the  estimated  fields  back  to  the  real  fields.  Nevertheless,  we  did  not  scale  measured 
pressure  and  currents  for  two  reasons.  Firstly,  non-uniformly  scaled  pressure  and 
current  data  will  no  longer  be  in  geostrophic  balance,  resulting  in  much  more  complex 
cross-covariance  functions  between  pressure  and  current.  Secondly,  in  contrast  to  the 
near-surface  Gulf  Stream  data,  the  deep  pressures  and  currents  have  only  small  spatial 
changes  in  eddy  variance;  thus  scaling  by  the  variance  should  not  be  necessary. 

•  isotropy.  Anisotropy  of  the  correlation  functions  in  this  region  is  not  justified  from 
our  previous  experience  with  mapping  the  thermocline  depth  fields  [8]. 

•  stationarity.  From  the  time  series  plots  of  measured  pressures  and  currents  (see  Fig¬ 
ure  9,  10  and  11  in  Appendix  C),  we  see  that  overall  the  time  series  are  approximately 
stationary. 

The  assumptions  above  have  reduced  the  problem  of  estimating  covariance  functions  to 
one  of  estimating  the  correlation  length  of  the  pressure  field.  Fitting  a  Gaussian  function  to 
data  covariance  estimates  in  Figure  3  results  in  a  correlation  length  of  about  90  kilometers. 
In  doing  the  least  square  fitting,  the  contribution  of  each  pair  of  measurement  sites  is 
weighted  by  the  length  of  the  concomitant  pressure  records  on  both  sites;  the  noise  variance 
in  the  data  cross-covariance  is  assumed  to  be  10  percent.  Also  shown  in  Figure  3  is  the  bin 
averages  of  the  estimated  data  cross-covariance.  Only  data  cross-covariance  estimated  by 
two  pressure  records  longer  than  30  days  are  used  in  bin  averages.  The  length  of  each  error 
bar  is  twice  as  big  as  the  standard  deviation  of  the  bin  average. 
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Distance  (km) 

Figure  3:  Spatial  correlation  function  of  the  pressure  fields. 

In  addition,  the  choice  for  the  values  of  the  noise  variance  and  correlation  length  in 
Table  6  were  also  based  on  the  following  experiment.  We  objectively  estimated  the  currents 
on  each  site,  first  using  ail  available  current  measurements  and  then  using  only  measure¬ 
ments  other  than  the  site  chosen.  In  both  cases,  all  combinations  of  noise  variance  (0.05, 
0.10)  and  correlation  lengths  from  50  km  through  130  km  (at  10  km  interval)  were  tried  in 
the  OA  procedure.  The  estimated  currents  were  compared  with  the  measured  ones  and  the 
resulting  root  mean  square  errors  were  compared  with  OA  predicted  errors.  The  parameter 
values  listed  in  Table  6  produced  best  overall  performance. 

The  covariance  matrices  A  and  C  in  (14)  and  (15)  are  calculated  from  (18)  and  (19). 
In  the  last  two  equations,  the  estimated  signal  covariance  functions  are  used  to  calculate 
the  elements  of  matrices  Aq  and  Co;  the  noise  variance  estimates  are  used  for  the  diagonal 
elements  of  E. 

3.2  Estimating  the  mean  fields 

Both  measured  and  estimated  fields  in  OA  formulae  (14)  —  (17)  are  assumed  to  have 
zero  mean.  Accordingly,  the  measured  field  should  be  demeaned  before  applying  the  OA 
procedure,  and  the  estimated  perturbation  field  should  have  mean  fields  restored  afterwards. 
This  subsection  describes  how  we  determine  those  mean  fields. 
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In  essence,  we  estimate  the  mean  pressure  fields  from  temporally  averaged  current  mea¬ 
surements  by  objectively  mapping  the  stream  function.  Thus  we  need  to  know  the  noise 
variance  in  the  averaged  measurements  and  the  covariance  functions  of  the  average  fields. 
In  this  case,  our  best  choice  is  that  the  noise  variance  is  five  percent  of  the  signal  variance 
for  both  averaged  pressure  and  averaged  currents  and  the  correlation  length  of  the  averaged 
pressure  fields  is  150  kilometers  (Table  6).  The  choice  was  guided  by  the  experiment  that 
we  mentioned  near  the  end  of  section  3.1,  in  which  we  found  that  our  choice  is  ‘best’  in  the 
sense  that  it  minimizes  the  difference  between  the  mean  estimated  currents  and  the  mean 
measured  ones. 

3.3  Data  flow  diagrams 

The  following  subsections  describe  the  processing  procedures  that  were  applied  to  the  mea¬ 
sured  pressures  and  currents  to  obtain  a  time  series  of  pressure  maps.  Figure  4  summarizes 
the  processing  steps. 

3.3.1  Dedrifting 

Several  bottom  residual  pressure  records,  H3  in  year  1,  and  G2,  H2,  H3,  II  in  year  2, 
appeared  to  have  trends  of  about  0. 1-0.2  decibar  in  one  year.  These  records  were  dedrifted 
by  fitting  an  exponential-plus-constant  curves  or  straight  lines  to  the  residual  pressure 
records  by  least  squares  method  and  subtracting  these  curves  from  the  records.  In  addition, 
there  was  a  slight  offset  between  the  two  G2  records  during  the  overlapping  period  of  the 
two  deployment  years.  This  offset  was  removed,  at  least  partly,  by  adding  half  the  average 
difference  (about  0.01  dbar)  to  the  record  for  year  2  and  subtracting  the  same  amount  from 
the  record  for  year  1.  Plots  in  Figure  5  and  6  show  all  the  relevant  PIES  records  before  and 
after  dedrifting. 

3.3.2  Reference  pressure  levels 

In  order  to  use  residual  pressure  records  from  PIES  pressure  measurements  in  making  the 
absolute  pressure  maps  (relative  to  an  overall  constant  independent  of  time  and  space), 
we  need  to  determine  the  reference  levels  for  the  residual  pressure  records.  Although  the 
residual  pressure  records  themselves  are  resolved  to  0.001  decibars,  the  directly  measured 
reference  levels  for  these  records  are  only  accurate  to  1.5  decibar  (page  11  in  [3])  due  to  the 
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Figure  4:  Data  processing  flowchart 
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Figure  6:  Dedrifted  PIES  records. 
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absolute  calibration  accuracy.  This  uncertainty  is  much  larger  than  the  standard  deviations 
of  these  pressure  records,  which  are  less  than  0.1  decibars.  Therefore,  we  develop  in  this 
section  a  better  way  of  determining  these  reference  levels. 

In  essence,  the  mean  pressure  held  mapped  from  current  records  was  used  as  reference 
level  for  the  PIES  measured  pressure  records.  The  following  scheme  was  used  to  ensure 
that  the  mean  reference  level  calculated  for  a  site  has  the  correct  average  period  consis¬ 
tent  with  measurement  period  of  that  site.  In  the  first  step,  diagrammed  at  the  top  of 
Figure  4,  objective  analysis  was  used  on  the  current  measurements  alone  to  map  pressure 
fields  ("u2p’’).  In  particular,  they  were  interpolated  to  the  PIES  instrument  sites.  Note, 
however,  that  while  the  pressure  fields  objectively  mapped  from  (u,v)  alone  are  spatially 
consistent,  the  reference  pressures  vary  temporally  from  map  to  map.  On  the  other  hand, 
the  bottom  pressure  record  each  has  a  temporally  consistent  reference  but  the  reference  can 
vary  from  site  to  site.  Let  P{t,s)  and  Po(f)  denote  the  interpolated  pressure  time  series  at 
site  s  from  current  meters  and  its  unknown  offset  from  the  absolute  pressure.  Let  p(t,s) 
and  po(s)  represent  the  measured  residual  pressure  from  PIES  at  site  s  and  its  unknown 
reference  level.  Since  the  absolute  pressure  at  a  given  time  and  site  is  unique,  we  have 

p{t,s)  +  Po{s)  -  P{t,  s)  +  Po{t). 

The  unknown  Po{t)  may  be  eliminated  by  subtracting  the  above  equation  for  any  one  site 
(we  used  G2)  from  the  equations  for  all  other  sites,  giving 

A,po(«)  =  A,P(f,s)  -  A,p((,s), 

where  is  the  spatial  difference  operator  with  respect  to  G2;  for  example,  A,p(t,  s)  ^ 
p(t,s)-p(t,G2). 

The  left  hand  term,  A«po(s),  of  the  above  equation  should  not  be  a  function  of  time. 
However,  the  right  hand  side,  estimated  from  real  data,  varies  with  time  due  to  errors 
(principally  those  in  mapping  P{t,s)  each  day  and  from  vertical  shear  changes  in  p(f,s) 
above  each  PIES).  In  the  light  of  this,  the  reference  level  for  each  site  s  on  the  left  hand 
side  of  the  above  equation  was  calculated  by  the  time  average  of  the  right  hand  side. 

3.3.S  Multivariate  mapping 

The  next  step,  midway  down  Figure  4,  is  to  combine  the  absolute  pressure  records  (Pb) 
and  current  measurement  records  (u,v)  in  mutivariate  objective  analysis,  to  obtain  pres- 
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sure  maps.  With  (Pb,  u,v)  inputs,  mapping  errors  are  reduced  because  pressure  and  its 
gradient  are  both  specified.  We  designate  these  maps  ‘‘preliminary”  because  of  an  iterative 
refinement  that  we  describe  next. 

The  objective  analysis  procedure  requires  that  the  input  data  have  the  mean  removed 
first.  Accordingly,  we  calculated  the  26  month  temporal  mean  pressure  field  and  the  corre¬ 
sponding  geostrophic  mean  current  field.  The  bottom  of  Figure  4  shows  (a)  we  subtracted 
these  mean  fields  from  the  absolute  pressure  records  and  the  current  records  to  get  “de¬ 
meaned”  records,  (b)  re-ran  the  multivariate  OA  program  on  the  demeaned  records,  and 
(c)  restored  the  mean  pressure  field  to  obtain  final  absolute  pressure  maps. 


A 
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4  Evaluation  of  the  maps 

This  section  assesses  the  precision  of  the  objective  pressure  maps  obtained  in  section  3  by 
comparing  the  time  series  derived  from  the  maps  and  their  statistics  with  those  from  the 
measurements. 

4.1  Comparison  of  mapped  fields  with  measured  fields 

In  this  section,  the  objectively  estimated  pressure  time  series  at  the  PIES  sites  are  com¬ 
pared  with  the  “measured”  absolute  pressures  as  determined  in  section  3.3.2.  Likewise,  the 
objectively  estimated  current  time  series  at  the  current  meter  sites  are  compared  with  the 
current  measurements.  The  comparison  plots  are  shown  in  Appendix  C  and  the  comparison 
statistics  are  tabulated  in  Tables  7,  8  and  9. 

The  comparison  plots  in  Appendix  C  show  good  agreement  between  measured  and 
mapped  pressure  and  currents. 

In  Tables  7,  p,  Pmea$  and  poa  denote  respectively  the  true  (that  is,  error  free)  pressure, 
the  PIES  measured  pressure  and  objectively  mapped  pressure.  On  most  sites,  the  correlar 
tion  coefficients  pp  between  pmeat  and  poa  are  above  90  percent.  In  carrying  out  the  OA 
procedure,  we  assumed  that  the  std(p)  is  a  constant  independent  of  sites.  If  we  choose 
std(p)  to  be  a  typical  value  of  0.7  decibar,  the  OA  estimated  std(poa  -  p)  agree  well  with 
“measured”  rms(poa  -  Pmeat)-  The  magnitude  of  the  difference  poa  -  Pmeat  is  in  naost  cases 
less  than  10  percent  of  std(pmea«)  as  shown  in  the  bias  column  in  Tables  7.  Note  the  bias  is 
defined  here  as  instead  of  the  more  widely  used  fn-Pmu  because  some  Pmeat'^ 

are  effectively  zero  for  the  given  averaging  periods. 

In  Table  8  and  9,  («,w),  and  (Uoa.Voo)  denote  respectively  the  true  (that 

is,  error  free)  current,  the  measured  current  and  objectively  mapped  current,  the  correlation 
coefficients  between  measured  and  mapped  currents  are  ail  90  percent  or  over.  It  appears, 
however,  that  the  OA  procedure  ‘smoothed’  out  some  variance  in  small  scale  motion  such 
that  when  8td(umeaj)  is  large  std(uoa)  tends  to  underestimate  it  and  tends  to 

underestimate 

4.2  Comparison  of  predicted  fields  with  measured  fields 

In  this  section,  we  examine  how  well  the  percentage  error  predicted  by  OA  procedure  reflects 
the  real  one.  To  do  this,  we  pick  one  pressure  measurement  site  at  a  time  and  remove  the 
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Table  7:  Comparison  between  measured  pressure  and  mapped  pressure.  (The  pressure  is  in 
decibar  unit  and  Len  denotes  record  length  in  days) _ 


Site 

Len 

Pmeat 

Poa 

std(p,„ea*) 

std(po„) 

bias 

Pp 

std(po.-p) 

std(j)) 

yig2 

472 

-0.01 

-0.01 

0.04 

0.04 

0.08 

0.93 

0.37 

0.23 

yig3 

365 

-0.02 

-0.02 

0.04 

0.04 

-0.08 

0.95 

0.30 

0.23 

ylh2 

369 

0.01 

0.01 

0.03 

0.03 

-0.01 

0.93 

0.37 

0.24 

ylh3 

427 

-0.02 

-0.02 

0.04 

0.03 

-0.06 

0.92 

0.40 

0.18 

ylh4 

424 

-0.04 

-0.04 

0.04 

0.04 

-0.13 

0.91 

0.41 

0.19 

ylhS 

430 

-0.08 

-0.08 

0.08 

0.06 

-0.04 

0.92 

0.41 

0.19 

ylhe 

374 

-0.08 

-0.07 

0.06 

0.05 

-0.05 

0.91 

0.42 

0.22 

ylil 

250 

0.00 

0.00 

0.04 

0.03 

0.05 

0.73 

0.69 

0.21 

yli2 

352 

-0.03 

-0.03 

0.04 

0.03 

-0.02 

0.84 

0.54 

0.22 

yli4 

370 

-0.08 

-0.08 

0.10 

0.08 

-0.04 

0.95 

0.35 

0.20 

yli5 

373 

-0.06 

-0.07 

0.09 

0.08 

0.02 

0.98 

0.24 

0.24 

y2g2 

438 

0.01 

0.01 

0.04 

0.04 

-0.05 

0.92 

0.38 

0.23 

y2g3 

445 

0.02 

0.02 

0.05 

0.05 

0.03 

0.98 

0.23 

0.24 

y2h2 

358 

0.04 

0.04 

0.03 

0.03 

0.02 

0.90 

0.46 

0.21 

y2h3 

362 

0.01 

0.01 

0.03 

0.03 

0.03 

0.87 

0.49 

0.18 

y2h4 

365 

0.00 

0.00 

0.06 

0.05 

0.04 

0.95 

0.30 

0.15 

y2h5 

366 

-0.01 

-0.02 

0.09 

0.09 

0.08 

0.99 

0.16 

0.18 

y2h6 

429 

-0.03 

-0.03 

0.08 

0.08 

0.01 

0.99 

0.17 

0.25 

y2il 

441 

0.03 

0.03 

0.04 

0.03 

0.05 

0.90 

0.43 

021 

y2i2 

438 

0.00 

0.00 

0.04 

0.03 

0.05 

0.87 

0.49 

0.18 

y2i3 

428 

-0.02 

-0.02 

0.08 

0.07 

0.00 

0.97 

0.26 

0.17 

y2i4 

424 

-0.04 

-0.04 

0.11 

0.10 

0.01 

0.99 

0.11 

0.17 

y2i5 

429 

-0.05 

-0.05 

0.10 

0.09 

0.00 

0.99 

0.17 

0.23 

Table  8:  Comparison  between  measured  currents  and  mapped  currents  (u  component).  (The 
current  is  in  meter /sec.  Len  denotes  record  length  in  days,  bias  is  defined  as 


denotes  correlation  coefficient  between  Umea»  ^d  Upg.) 


Site 

Len 

Umeo* 

Uoa 

Std(tt,„ea.) 

std(«o«) 

bias 

Pu 

rnis(uoa  “ttffWM ) 

yig2 

371 

0.01 

0.00 

0.06 

0.05 

0.10 

0.97 

0.26 

0.27 

yig3 

137 

-0.01 

-0.02 

0.05 

0.05 

0.14 

0.96 

0.28 

0.26 

ylh4 

354 

-0.06 

-0.05 

0.10 

0.08 

-0.04 

0.96 

0.32 

0.23 

ylhS 

370 

-0.03 

-0.04 

0.08 

0.06 

0.07 

0.90 

0.44 

0.23 

yli2 

352 

-0.06 

-0.06 

0.09 

0.08 

0.01 

0.96 

0.30 

0.26 

yli3 

146 

0.00 

-0.01 

0.09 

0.07 

0.09 

0.92 

0.40 

0.22 

yli4 

87 

0.04 

0.05 

0.04 

0.04 

-0.36 

0.93 

0.38 

0.23 

y2g2 

438 

0.02 

0.02 

0.09 

0.08 

0.01 

0.98 

0.23 

0.27 

y2h2 

805 

-0.07 

-0.07 

0.08 

0.06 

-0.01 

0.98 

0.24 

0.27 

y2h3 

437 

-0.03 

-0.03 

0.10 

0.08 

0.05 

0.94 

0.37 

0.23 

y2h5 

427 

-0.03 

-0.03 

0.11 

0.10 

-0.02 

0.97 

0.26 

0.22 

y2h6 

365 

0.01 

0.01 

0.06 

0.06 

0.00 

0.97 

0.23 

0.28 

y2il 

806 

-0.06 

-0.06 

0.08 

0.07 

0.00 

0.97 

0.26 

0.28 

y2i2 

440 

-0.05 

-0.05 

0.11 

0.09 

-0.01 

0.95 

0.33 

0.23 

y2i3 

348 

-0.04 

-0.03 

0.13 

0.14 

-0.05 

0.97 

0.24 

0.21 

y2i4 

353 

0.00 

0.00 

0.10 

0.09 

0.00 

0.94 

0.33 

0.22 

y2i5 

805 

0.04 

0.04 

0.10 

0.10 

0.01 

0.99 

0.16 

0.28 

y2ml3 

354 

-0.01 

-0.01 

0.13 

0.11 

0.02 

0.97 

0.26 

0.20 

Table  9:  Comparison  between  measured  currents  and  mapped  currents  (v  component).  (See 
caption  of  Table  8  for  notational  convention) _ 


Site 

Len 

^meat 

Voa 

Std(Vmea.) 

std(i;oa) 

bias 

Pv 

8td(voa-v) 

Std(v) 

yig2 

371 

0.00 

0.00 

0.06 

0.05 

-0.01 

0.97 

0.25 

0.27 

yig3 

137 

-0.04 

-0.06 

0.08 

0.07 

0.15 

0.97 

0.24 

0.26 

ylh4 

354 

-0.05 

-0.05 

0.08 

0.07 

-0.01 

0.96 

0.30 

0.23 

ylh5 

370 

-0.05 

-0.05 

0.10 

0.10 

0.05 

0.97 

0.26 

0.23 

yli2 

352 

-0.01 

-0.01 

0.05 

0.04 

-0.02 

0.95 

0.40 

0.26 

yli3 

146 

0.03 

0.04 

0.08 

0.08 

-0.12 

0.97 

0.25 

0.25 

yli4 

87 

0.11 

0.12 

0.09 

0.09 

-0.09 

0.98 

0.19 

0.25 

y2g2 

438 

0.01 

0.01 

0.08 

0.07 

-0.01 

0.98 

0.25 

0.28 

y2h2 

805 

-0.01 

-0.01 

0.05 

0.04 

-0.06 

0.93 

0.39 

0.26 

y2h3 

437 

-0.01 

-0.01 

0.07 

0.05 

-0.03 

0.91 

0.45 

0.21 

y2h5 

427 

-0.05 

-0.05 

0.11 

0.11 

0.00 

0.98 

0.20 

0.24 

y2h6 

365 

-0.06 

-0.05 

0.09 

0.09 

-0.04 

0.98 

0.18 

0.27 

y2il 

806 

-0.01 

-0.01 

0.05 

0.04 

-0.02 

0.91 

0.43 

0.27 

y2i2 

440 

-0.02 

-0.02 

0.06 

0.04 

0.00 

0.91 

0.42 

0.25 

y2i3 

348 

-0.05 

-0.04 

0.09 

0.08 

-0.06 

0.94 

0.33 

0.22 

y2i4 

353 

-0.03 

-0.03 

0.10 

0.09 

0.04 

0.99 

0.18 

0.26 

y2i5 

805 

0.04 

0.04 

0.12 

0.12 

0.00 

0.99 

0.13 

0.28 

y2ml3 

354 

•0.04 

-0.05 

0.10 

0.10 

0.08 

0.98 

0.22 

0.21 
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pressure  measurement  record  on  the  site  and  run  the  OA  procedure  as  if  that  pressure 
record  did  not  exist.  As  a  result,  we  obtain  an  interpolated  pressure  record  and  a  predicted 
percentage  error  on  the  measurement  site.  Then  we  compare  the  observed  pressure  record 
with  the  interpolated  one  and  calculate  an  observed  percentage  error.  The  results,  along 
with  other  comparison  statistics,  are  displayed  in  Table  10.  Similar  results  for  currents  are 
shown  in  Table  11  and  Table  12. 

It  may  be  worth  pointing  out  that  this  is  an  extreme  test  since  it  leaves  such  a  big  hole 
in  the  array.  Most  places  were  surrounded  by  closer  measurements  than  in  such  a  case.  The 
predicted  errors  (last  column)  in  Tables  10,  11  and  Table  12  tend  to  be  smaller  than  the 
measured  errors  (second  to  the  last  column).  Spedhcally,  the  discrepancy  seems  worst  on 
array  edges,  less  in  the  interior  and  least  at  site  M13.  So  the  fewer  nearby  measurements 
the  worse  discrepancy.  This  may  be  due  to  the  approximate  nature  of  the  assumptions 
in  section  3.1.  There  may  be,  however,  another  cause  for  the  discrepancy.  In  the  case  of 
pressure  estimate,  this  can  be  seen  as  foUows.  In  estimating  Poa  on  a  site,  Pmea*  on  the 
site  was  taken  out.  Thus  the  measurement  error  in  pmea»  should  be  uncorrelated  to  the 
mapping  error  in  poa-  Accordingly,  we  have 

rms(Poa  -  Pmeat )  rms(poa  -  p  +  p  -  Pmea.)  _  rms(poa  “  p)  .  rm8(p  -  Pme«.) 

Std(pniea<)  Std(pnieo»)  Std(pnieo»)  8^*^(Pmeo#) 

The  first  and  second  terms  on  the  right  can  measure  respectively  the  mapping  and  measure¬ 
ment  errors.  Therefore,  overestimate  the  mapping  error  in  po*  given  by 

Std(poa-p) 

std(p)  ■ 

4.3  Checking  pressure  maps  against  measured  currents 

The  geostrophic  currents  from  pressure  maps  are  compared  with  current  measurements  in 
section  5.2  by  overlaying  the  observed  current  vectors  on  the  daily  pressure  maps.  Like¬ 
wise,  the  mean  geostrophic  currents  from  pressure  maps  are  compared  with  mean  current 
measurements  in  section  5.1  by  overlaying  the  mean  observed  current  vectors  on  the  mean 
pressure  maps.  From  visiual  inspection,  the  current  vectors  in  these  plots  almost  always 
run  parallel  with  the  pressure  contour  lines. 
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Table  10:  Comparison  between  measured  pressure  and  predicted  pressure.  (See  caption  of 
Table  7  for  notational  convention) _ _ 


Site 

Len 

Pmeaa 

Poa 

8t<i(Pmettj) 

std(poo) 

bias 

Pp 

rms(i>M-Dm«u) 

Std(pm«a«) 

8td(p«i-p) 

8td(j)) 

yig2 

472 

-0.01 

-0.01 

0.04 

0.04 

-0.13 

0.71 

0.78 

0.34 

yig3 

365 

-0.02 

-0.01 

0.04 

0.04 

0.15 

0.78 

0.64 

0.34 

ylh2 

369 

0.01 

0.01 

0.03 

0.03 

0.04 

0.59 

0.92 

0.38 

ylh3 

427 

-0.02 

-0.02 

0.04 

0.03 

0.08 

0.81 

0.59 

0.22 

ylli4 

424 

-0.04 

-0.03 

0.04 

0.04 

0.20 

0.79 

0.63 

0.24 

ylhS 

430 

-0.08 

-0.08 

0.08 

0.06 

0.06 

0.78 

0.63 

0.26 

ylh6 

374 

-0.08 

-0.07 

0.06 

0.06 

0.09 

0.58 

0.87 

0.31 

ylil 

250 

0.00 

0.00 

0.04 

0.04 

-0.11 

0.17 

1.24 

0.30 

yli2 

352 

-0.03 

-0.03 

0.04 

0.04 

0.08 

0.47 

0.96 

0.31 

yli4 

370 

-0.08 

-0.08 

0.10 

0.07 

0.08 

0.84 

0.55 

0.27 

yliS 

373 

-0.06 

-0.06 

0.09 

0.06 

0.01 

0.86 

0.54 

0.36 

y2g2 

438 

0.01 

0.01 

0.04 

0.04 

0.13 

0.64 

0.89 

0.35 

y2g3 

445 

0.02 

0.01 

0.05 

0.05 

-0.09 

0.83 

0.56 

0.39 

y2h2 

358 

0.04 

0.04 

0.03 

0.04 

-0.04 

0.74 

0.82 

0.28 

y2h3 

362 

0.01 

0.01 

0.03 

0.03 

-0.05 

0.71 

0.75 

0.23 

y2h4 

365 

0.00 

0.00 

0.06 

0.05 

-0.05 

0.92 

0.40 

0.18 

y2h5 

366 

-0.01 

-0.02 

0.09 

0.09 

-0.11 

0.96 

0.27 

0.23 

y2h6 

429 

-0.03 

-0.04 

0.08 

0.07 

-0.04 

0.88 

0.47 

0.42 

y2il 

441 

0.03 

0.03 

0.04 

0.04 

o 

o 

1 

0.67 

0.80 

0.28 

y2i2 

438 

0.00 

0.00 

0.04 

0.04 

-0.07 

0.72 

0.73 

0.23 

y2i3 

428 

-0.02 

-0.02 

0.08 

0.07 

-0.02 

0.93 

0.38 

0.20 

y2i4 

424 

-0.04 

-0.04 

0.11 

0.10 

-0.02 

0.99 

0.17 

0.21 

y2i5 

429 

-0.05 

-0.05 

0.10 

0.09 

-0.01 

0.94 

0.35 

0.34 
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Table  11:  Comparison  between  measured  currents  and  predicted  currents  (u  component). 
(See  caption  of  Table  8  for  notational  convention) _ 


Site 

Len 

^meat 

Uoa 

Std(u„ea*) 

Std(Uoa) 

bias 

Pu 

rmS(uaa-Um<a<) 

8td(tt«n«aj) 

St(l(Uo«‘-1A) 

8td(u) 

yig2 

371 

0.01 

-0.01 

0.06 

0.05 

-0.35 

0.52 

0.90 

0.51 

yig3 

137 

-0.01 

-0.03 

0.05 

0.05 

-0.31 

0.58 

0.87 

0.49 

ylh4 

354 

-0.06 

-0.05 

0.10 

0.07 

0.10 

0.76 

0.65 

0.32 

ylhS 

370 

-0.03 

-0.04 

0.08 

0.07 

-0.09 

0.46 

0.98 

0.34 

yli2 

352 

-0.06 

-0.06 

0.09 

0.07 

0.02 

0.60 

0.81 

0.48 

yli3 

146 

0.00 

-0.01 

0.09 

0.07 

-0.13 

0.63 

0.80 

0.31 

yli4 

87 

0.04 

0.07 

0.04 

0.04 

0.76 

0.63 

0.84 

0.35 

y2g2 

438 

0.02 

0.01 

0.09 

0.07 

-0.04 

0.66 

0.76 

0.52 

y2h2 

805 

-0.07 

-0.06 

0.08 

0.05 

0.09 

0.48 

0.91 

0.56 

y2h3 

437 

-0.03 

-0.04 

0.10 

0.08 

-0.13 

0.63 

0.78 

0.34 

y2h5 

427 

-0.03 

-0.03 

0.11 

0.09 

0.07 

0.85 

0.53 

0.33 

y2h6 

365 

0.01 

0.00 

0.06 

0.08 

-0.08 

0.68 

0.98 

0.58 

y2il 

806 

-0.06 

-0.07 

0.08 

0.06 

-0.04 

0.28 

1.08 

0.59 

y2i2 

440 

-0.05 

-0.05 

0.11 

0.09 

0.04 

0.72 

0.71 

0.34 

y2i3 

348 

-0.04 

-0.03 

0.13 

0.15 

0.08 

0.92 

0.44 

0.28 

y2i4 

353 

0.00 

0.00 

0.10 

0.09 

0.00 

0.77 

0.66 

0.32 

y2i5 

805 

0.04 

0.03 

0.10 

0.10 

-0.04 

0.76 

0.68 

0.59 

y2ml3 

354 

-0.01 

o 

o 

0.13 

0.11 

-0.03 

0.90 

0.43 

0.26 

Table  12:  Comparison  between  measured  currents  and  predicted  currents  (v  component). 
(See  caption  of  Table  8  for  notational  convention) _ 


Site 

Len 

Vmeo* 

Voa 

std(t;„,eo.) 

std(roa) 

bias 

Pv 

rms(vM-vnMu) 

Std(VllMu) 

yig2 

371 

0.00 

0.00 

0.06 

0.06 

-0.01 

0.47 

1.02 

0.56 

yig3 

137 

-0.04 

-0.08 

0.08 

0.07 

-0.46 

0.65 

0.80 

0.51 

ylh4 

354 

-0.05 

-0.05 

0.08 

0.06 

0.04 

0.78 

0.63 

0.34 

ylh5 

370 

-0.05 

-0.06 

0.10 

0.12 

-0.15 

0.87 

0.58 

0.34 

yli2 

352 

-0.01 

-0.01 

0.05 

0.04 

0.07 

0.12 

1.13 

0.47 

yli3 

146 

0.03 

0.06 

0.08 

0.09 

0.33 

0.86 

0.61 

0.40 

yli4 

87 

0.11 

0.13 

0.09 

0.09 

0.23 

0.87 

0.50 

0.41 

y2g2 

438 

0.01 

0.01 

0.08 

0.07 

0.04 

0.29 

1.10 

0.60 

y2h2 

805 

-0.01 

0.00 

0.05 

0.05 

0.22 

0.17 

1.24 

0.50 

y2h3 

437 

-0.01 

-0.01 

0.07 

0.05 

0.04 

0.55 

0.85 

0.29 

y2h5 

427 

-0.05 

-0.05 

0.11 

0.10 

0.03 

0.87 

0.49 

0.38 

y2h6 

365 

-0.06 

-0.04 

0.09 

0.11 

0.16 

0.84 

0.67 

0.54 

y2il 

806 

-0.01 

0.00 

0.05 

0.06 

0.10 

0.02 

1.56 

0.52 

y2i2 

44U 

-0.02 

-0.01 

0.06 

0.05 

0.03 

0.31 

1.11 

0.42 

y2i3 

348 

-0.05 

-0.04 

0.09 

0.08 

0.12 

0.76 

0.66 

0.31 

y2i4 

353 

-0.03 

-0.04 

0.10 

0.08 

-0.12 

0.86 

0.52 

0.43 

y2i5 

805 

0.04 

0.04 

0.12 

0.11 

0.01 

0.80 

0.60 

0.58 

y2ml3 

354 

-0.04 

-0.05 

0.10 

0.10 

•0.14 

0.93 

0.39 

0.27 
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5  Bottom  Pressure  Maps 

This  section  displays  the  mean  and  daily  pressure  maps. 

5.1  Mean  pressure  maps 

Figure  7  shows  mean  bottom  pressure  maps  with  measured  (solid  line  vectors)  and  mapped 
(dash  line  vectors)  mean  currents.  Contour  levels  are  in  decibars.  The  scale  of  the  currents 
is  3  kilometer  in  the  maps  corresponding  to  1  centimeter  per  second.  The  mean  is  taken  over 
the  period  from  June  16,  1988  through  May  24,  1989  for  year  1,  and  from  August  29,  1989 
through  August  7, 1990  for  year  2.  Both  pressure  and  current  are  averaged  over  same  period. 
Sites  with  measured  pressure  and  current  records  contributing  to  the  maps  are  respectively 
indicated  by  o’s  and  x’s.  However,  mean  current  vectors  are  not  shown  on  three  current 
measurement  sites  in  year  1  because  the  measurements  spanned  only  part  of  the  averaging 
period  shown.  A  description  of  the  coordinate  system  is  contained  in  section  5.2. 
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5.2  Daily  pressure  maps  overlaid  with  current  measurements 

Bottom  pressure  maps  from  June  13,  1988  through  August  7,  1990  are  shown  at  daily 
interval  in  this  section.  The  original  maps  were  produced  at  six  hour  interval  but  only  the 
ones  at  1200  UT  are  displayed. 

The  overall  format  of  the  pressure  maps  is  designed  to  facilitate  comparison  with  the 
thermocline  depth  maps  in  the  Central  Array  for  the  same  time  period  [8].  To  this  end,  the 
border  of  the  thermocline  depth  maps  in  [8}  is  shown  here  with  each  pressure  map.  The 
region  of  the  pressure  maps  corresponds  to  the  boxed  area  in  Figure  1.  The  latitude  and 
longitude  of  the  four  corners  of  the  box  are  listed  in  Table  1.  The  coordinate  system  used 
here  is  the  same  as  that  used  for  the  thermocline  depth  maps  in  [8].  The  origin  is  located 
at  38**N,  68**W.  The  scaling  factors  converting  latitude  and  longitude  to  kilometers  are  111 
km  per  degree  latitude  and  82.263  km  per  degree  longitude.  The  grid  is  then  rotated  15” 
so  that  the  horizontal  axis  is  directed  075”T  and  vertical  axis  345”T. 

Shown  in  each  map  are  the  bottom  pressure  contours,  percentage  error  contours,  ob¬ 
served  current  vectors,  current  meter  mooring  positions  and  PIES  sites.  The  pressure  field 
is  contoured  at  0.015  decibar  intervals.  The  negative  and  positive  contour  lines  are  denoted 
respectively  by  dashed  and  solid  lines  to  tell  troughs  from  crests.  The  estimated  errors 
are  represented  as  percentage  of  the  standard  deviation  of  the  pressure  field,  which  is  0.6 
decibar  on  average.  Errors  are  contoured  at  3%  interval  with  those  over  48%  shown  by  solid 
contours  and  those  between  30%  and  48%  by  dashed  contours.  Errors  lower  than  30%  are 
not  contoured.  The  current  vectors  are  shown  by  solid  lines  with  3  kilometer  in  the  maps 
corresponding  to  1  centimeter  per  second.  Small  squares  and  marks  denote  respectively 
current  meter  positions  and  PIES  sites. 

Figure  8:  Daily  pressure  maps  overlaid  with  current  measurements  (pages  33  —  163).  See 
paragraphs  above  for  description. 
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Appendices 

A  OA  Program  List 

A.l  to  1*  OA  program 

function  Cpi,ai]«oall(xi.yi,xm.ym.pm,Vn.CLENGTH) 

Xf unction  Cpi>«i]*o%li(xi,yi.xm.ym.pm,Vn.CLENGTH) 

X  scalar  objsctivs  analysis 

X  pi  :  intsrpolatad  values  vith  columns  being  the  time  and 
X  rows  being  the  interpolation  points  (NI  by  NT) ; 

X  ei  :  estimated  error  in  percentage  of  signal  variance  (NI  by  1) ; 

X  xi,  yi:  x-y  coordinates  of  interpolation  points  (NI  by  1): 

X  xm,  ya:  x-y  coordinates  of  measurement  sites  (NM  by  1); 

X  pm  :  array  of  measured  values  vith  columns  being  the  time 

X  and  rows  being  the  measurement  sites  (NM  by  NT) : 

X  Vn  :  noise  variance  in  percentage  of  signal  variance  (scalar  or  1  by  NM) ; 
X  CLENGTH:  correlation  length  (scalar); 

X  units:  xi,yi,xm,ym  and  CLENGTH  should  have  the  same  units; 

X  input  pm  and  output  pi  should  have  the  same  units. 

X 

NM  •  length(xm); 

'  NI  «  length(xi) ; 

X  construct  correlation  matrices 

r  ■  sqrt(  (  xm*ones(l,NM)  -  ones(NN,l)*xm*  ).*2  ♦... 

(  ym*one8(l,NN)  -  one8(NM,l)*ym*  ).*2  ); 

A  ■  CFUNpp(r, CLENGTH) ;  X  pick  correlation  function  here. 

X  CFUNpp(r,L)  ■  exp(-(r/L) .“2)  ; 

A(1:NM*1:NM“2)»A(1:NM+1:NM*2) ,*(1  ♦  Vn);  X  taking  care  of  meaeurement  errora. 

r  ■  8qrt(  (  xi*one8(l,NN)  -  one8(NI,l)*xm*  ).“2  ♦... 

(  yisonead.NM)  -  one8(NI,l)*yn’  ).“2  ); 

C  ■  CFUMpp(r. CLENGTH) ;  X  pick  correlation  function  here. 

X  mapping 
pi  ■  Cs(A\pm); 

ei  ■  1  -  diag( (C/A) *C’)/CFUNpp(0, CLENGTH);  X  pick  correlation  function  here. 


Xnufflber  of  meaeurement  aitea. 
Xnufflber  of  interpolation  pointa. 
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A.2  *2  to  1’  OA  program 

function  Cpi.«i]*oa2lCxi,yi,zm.ym.uiB.vm,Vn,CLENGTH) 

Xf unction  Cpi.«i]«oa21(zi,7i.xiii.ym.uffl,vm,Vn,CLENGTH) 

X  "2  to  1"  objactiva  analysis 

X  pi  :  intarpolatad  valuas  with  columns  baing  tha  tima  and 
X  rows  baing  tha  intarpolation  points  (NI  by  NT) ; 

X  ai  :  astimatad  arror  in  parcantaga  of  signal  varianca  (NI  by  1) ; 

X  zi.  yi:  z-y  coordinatas  of  intarpolation  points  (NI  by  1): 

X  zm,  ym:  z-y  coordinatas  of  maasuramant  sitas  (NN  by  1) ; 

X  urn,  vm:  array  of  maasurad  valuas  with  columns  baing  tha  tima 
X  and  rows  baing  tha  maasuramant  sitas  (NM  by  NT) ; 

X  Vn  :  noisa  varianca  in  parcantaga  of  signal  varianca  (scalar  or  1  by  2*NM); 
X  CLENGTH:  corralation  langth  (scalar): 

X  units:  zi,yi,zm,ym  and  CLENGTH  should  hava  tha  sama  units; 

X  units  of  inputs  um.vm  and  output  pi  should  hava  CONSISTENT  units 

X  with  aach  othar  and  with  corralation  functions  CFUN??. 

NM  ■  langth(zm) ;  Xnumbar  of  currant  moorings . 

NI  *  langth(zi) :  Xnumbar  of  intarpolation  points. 

X  construct  corralation  matricas 

dz  ■  -  zm*onaa(l,NM)  ♦  onas(NM.l)*zm’ ;  Xdz(r,c)azm.c  -  zm.r; 
dy  ■  -  ymvonasd.NM)  ♦  onas(NM,l)*yB' ; 

Auu  a  CFUNuu (dz,dy, CLENGTH ) ;  Xpick  corralation  fun  hara. 

Avv  ■  CFUNvv (dz.dy, CLENGTH) ;  Xpick  corralation  fun  hara. 

Auv  ■  CFUNuv(dz, dy, CLENGTH) ;  Xpick  corralation  fun  hara. 

A*  CAuu ,  Auv ;  Auv ’ .  Avv] ; 

L  >  2*NH; 

A(1:L+1:L‘2)  »  A(1:L+1:L‘2)  .*  (1  ♦  Vn) ;  Xtaka  cara  of  maasuramant  arrors. 

dz  ■  -  zi*ona8(l,NM)  ♦  onas(NI,i)*zm’ ;  Xdz(r,c)»zm_c  -  zi.r 
dy  ■  -  yi«onas(l,NM)  ♦  onas(NI,l)*ym' ; 

Cpu  ■  CFUNpu(dz,dy, CLENGTH);  Xpick  corralation  fun  hara. 

Cpv  •  CFUHpv(dz,dy, CLENGTH);  Xpick  corralation  fun  hara. 

CaCCpu,  Cpv] : 

X  mapping 

pi  «  C«(A\[um:vm]) ; 

ai  ■  1  -  diag( (C/A) •C’)/CFUNpp(0, CLENGTH);  Xpick  corralation  fun  hara. 
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A.3  *2  to  2’  OA  program 

function  [ui,vi.«u.«T]«oa22(xi.yi,za.7m.ua.vB.Vn.CLENGTH} 

Xfunctlon  Cui.vi.«u.«v]>oa22(xl.yi.za.yB,uB,vm.Vn,CLENGTH) 

X  "2  to  2"  objoctivo  analysis 

X  ui.Ti  :  intarpolatsd  valuss  sith  colusms  bsing  ths  tins  and 
X  rovs  bsing  ths  intsrpolation  points  (NI  by  NT) ; 

X  su.sT  :  sstinatsd  srror  in  psrcsntags  of  signal  variancs  (VI  by  1); 

X  xi,  yi:  x-y  coordinatss  of  intsrpolation  points  (MI  by  1); 

X  XB,  ym:  x-y  coordinatss  of  nsasursnsnt  sitss  (NM  by  1) ; 

X  UB,  vb:  array  of  Bsasursd  valuss  sith  colunns  bsixig  ths  tiBS 
X  and  rovs  bsing  ths  BsasursBsnt  sitss  (NM  by  NT) ; 

X  Vn  :  noiss  variants  in  psrcsntags  of  signal  variants  (scalar  or  1  by  2*MM) ; 
X  CLENGTH:  corrslation  Isngth  (scalar); 

X  units:  xi.yi.xB.yB  and  CLENGTH  should  havs  ths  sans  units; 

X  units  of  inputs  ub.vb  and  output  ui.vi  should  havs  CONSISTENT  units 

X  vith  sach  othsr  and  vith  corrslation  functions  CFUN??. 

NM  >  Isngth(xB);  Xnuabsr  of  currant  Boorings. 

VI  ■  Isngth(xi);  Xnuabsr  of  intsrpolation  points. 

X  construct  corrslation  Batricss 

dx  ■  -  XBSonssd.NM)  ♦  onss(NH,l)*ZB' ;  X  dx(r,c)  •  xb.c  -  XB_r 
dy  ■  -  yB*onss(l.NM)  *  onss(NM,l)*yiB* ; 

Auu  «  CFUNuu (dx.dy, CLENGTH ) ;  Xpick 

Aw  «  CFUVvv (dx.dy, CLENGTH);  Xpick 

Auv  >  CFUNuv(dx.dy, CLENGTH);  Xpick 

A" [Auu,  Auv;  Auv*,  Avv] ; 

L  •  2*HN; 

A(1:L+1:L“2)  ■  A(1:L+1:L‘2)  .*  (1  ♦  Vn); 

dx  ■  -  xi*onss(l,NM)  *  onss(NI,l)*XB’ ; 
dy  ■  -  yi*onss(l.NH)  ♦  onss(NI,l)*yB* ; 

Cuu  »  CFUNuu(dx.dy. CLENGTH);  Xpick 

Cuv  ■  CFUVuv (dx.dy, CLENGTH ) ;  Xpick 

Cvu  -  CFUMvu(dx.dy, CLENGTH):  Xpick 

Cw  ■  CFUNvv(dx,dy, CLENGTH);  Xpick 

CuaCCuu,  Cuv]; 

Cv« [Cvu ,  Cvv] ; 

X  Bapping 

ui  ■  Cu*(A\[ub;vb]) : 
vi  •  Cvs(A\Cub;vb]); 

stt  ■  1  -  diaf((Cu/A)*Cu')/CFUHuu(0.0. CLENGTH);  Xpick  corrslation  fun  hsrs. 
sv  -  I  -  diag((Cv/A)sCv*)/CFUNvv(0,0,CLENGTH);  Xpick  corrslation  fun  hsrs. 


corrslation  fun  hsrs. 
corrslation  fun  hsrs. 
corrslation  fun  hors. 


Xtaks  cars  of  Bsasuroasnt  orrors. 

Xdx(r,c)«XB_c  -  xi_r 

corrslation  fun  hsrs. 
corrslation  fun  hsrs. 
corrslation  fun  hsrs. 
corrslation  fun  hsrs. 
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A.4  ‘3  to  1*  OA  program 

function  Cpi,®i]»oa31(xi.yi,xp,yp.p,xv,yv,u,v,Vn,CLENGTH) 

Xf unction  Cpi.ei]«oa31(xl,yi,xp.yp.p,xv.yv.u,T,Vn,CLENGTH) 

X  "3  to  1“  objective  analysis 

X  pi  :  interpolated  values  vith  columns  being  the  time  and 
X  roes  being  the  interpolation  points  (HI  by  NT) : 

X  ei  :  estimated  error  in  percentage  of  signal  variance  (HI  by  1) ; 

X  xi,  yi:  x-y  coordinates  of  interpolation  points  (NI  by  1); 

X  xp,  yp:  x-y  coordinates  of  pressure  measurement  sites  (NP  by  1); 

X  p  :  measured  pressures  eith  columns  being  the  time 

X  and  rous  being  the  measurement  sites  (NP  by  NT) : 

X  XV.  yv:  x-y  coordinates  of  current  measurement  sites  (NV  by  1); 

X  u,  V  :  measured  currents  eith  columns  being  the  time 
X  and  roes  being  the  measurement  sites  (NV  by  NT) ; 

X  Vn  :  noise  variance  in  percentage  of  signal  variance 
X  (scalar  or  1  by  NP+2*NV); 

X  CLENGTH:  correlation  length  (scalu); 

X  units:  xi.yi.xp.yp.xv.yv  and  CLENGTH  should  have  the  same  units; 

X  units  of  inputs  p,u,v  and  output  pi  should  be  CONSISTENT 

X  eith  each  other  and  eith  the  correlation  functions  CFUN??. 

NP  ■  length(xp);  Xnumber  of  pressure  gauges. 

NV  ■  length (xv);  Xnumber  of  current  moorings. 

NI  •  length(xi);  Xnumber  of  interpolation  points. 

X  construct  correlation  matrices 

Rpp  ■  sqrt(  (  xp*ones(l,NP)  -  ones(NP,l)*xp’  ).‘2  ♦ — 

(  yp*ones(l,NP)  -  ones(NP,l)*yp’  ).‘2  ); 
dxpv  ■  -  xp*ones(l,NV)  ♦  ones(NP,l)*xv» ;  X  dxpv(r,c)  •  xv.c  -  xp.r 
dypv  ■  -  yp*ones(l,NV)  +  one8(NP,l)*yv’ ; 
dxvv  «  -  xv*ones(l,NV)  *  ones(NV,l)*xv’ ; 
dyvv  ■  -  yv*ones(l  ,NV)  ♦  ones(NV,l)'»yv' ; 

App  •  CFUMpp (Rpp, CLENGTH) ; 

Apu  •  CFUNpu( dxpv, dypv, CLENGTH) 

Apv  ■  CFUNpv(dxpv, dypv, CLENGTH) 

Auu  ■  CFUNuu(dxvv, dyvv, CLENGTH) 

Avv  ■  CFUNvv (dxvv, dyvv, CLENGTH) 

Auv  ■  CFUNuv( dxvv, dyvv, CLENGTH) 

A«C  App,  Apu,  Apv; . . . 

Apu', Auu,  Auv;... 

Apv’ ,Auv’ ,Avv] ; 

L  *  NP«2*MV; 

A(1:L+1:L“2)  ■  A(1 :L*1 :L“2) .•(!  ♦  Vn);  Xtaking  care  of  measurement  errors. 


X  pick  correlation  fun  here. 
X  pick  correlation  fun  here. 
X  pick  correlation  fun  here. 
X  pick  correlation  fun  here. 
X  pick  correlation  fun  here. 
X  pick  correlation  fun  here. 
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Rip  ■  «qrt(  (  zi*on«8(l,MP}  -  on«8(NI.l)*zp*  ),'2  ♦... 

(  yi*on«8(l,IIP)  -  on«8(NI.l)*yp»  ),“2  ); 
dziv  •  -  zi*on«8(l.llV)  ♦  on«8(NI,l)*zv* ;  Xdx(r,c)«zv.c  -  zi_r 

dyiv  ■  -  yl*on«8(l,NV)  ♦  on«8(NI.l)*yT* ; 

Cpp  ■  CFUNppCRip.CLENGTH) :  X  pick  corralation  fun  h«r«. 

Cpu  ■  CFUNpu(dziv,dyiT,CLENQTH);  X  pick  corr«lation  fun  h«r8. 

Cpv  ■  CFUHpvCdziv.dyiv.CLEHGTH) ;  X  pick  correlation  fun  here. 

C«CCpp,  Cpu,  Cpv]; 

X  mapping 

pi  ■  C*(A\[p;u;v]) : 

ei  ■  1  -  diag((C/A)*C’ )/CFUMpp(0,CLENGTH) ;  Xpick  correlation  fun  hare. 
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A.5  ‘3  to  2’  OA  program 

function  Cui,vi,«u,«v]«oa32(xi,yi,xp,yp,p,xv,yv,u,v,Vn,CLENGTH) 

X  Cui,vi,*u,«v]«oa32(xi.yi,xp,yp.p,xv,yv,u.T,Vn,CLEHGTH) 

X  "3  to  2"  objoctiva  analysis 

X  ui.vi  :  intarpolatad  valuas  aith  coluoms  baing  tha  tima  and 
X  rows  baing  tha  intarpolation  points  (HI  by  HT) ; 

X  au.ar  :  astioatad  error  in  parcantaga  of  signal  varianca  (HI  by  1); 

X  xi,  yi:  x-y  coordinatas  of  intarpolation  points  (HI  by  1); 

X  xp,  yp:  x-y  coordinatas  of  prasstira  maasuramant  sitas  (HP  by  1); 

X  p  :  naasurad  prassuras  with  columns  baing  tha  tima 

X  and  rows  baing  tha  maasuramant  sitas  (HP  by  HT) ; 

X  xw.  yv:  x-y  coordinatas  of  currant  maasuramant  sitas  (HV  by  1); 

X  u,  V  :  maasurad  currants  with  columns  baing  tha  tima 
X  and  rows  baing  tha  maasuramant  sitas  (HV  by  HT) ; 

X  Vn  :  noisa  varianca  in  parcantaga  of  signal  varianca 
X  (scalar  or  1  by  HP-»2*HV) ; 

X  CLEHGTH:  corralation  langth  (scalar); 

X  units:  xi.yi.xp.yp.xv.yv  and  CLEHGTH  should  hava  tha  sama  units; 

X  units  of  inputs  p,u,v  and  output  ui.vi  should  ba  COHSISTEIT 

X  with  aach  othar  and  with  tha  corralation  functions  CFUH??. 

HP  «  langth(xp);  Xnumbar  of  prassuro  gaugas. 

HV  ■  langth(xv) ;  Xnumbar  of  currant  moorings. 

HI  ■  langth(xi) ;  Xnumbar  of  intarpolation  points. 

X  construct  corralation  matricas 

Rpp  ■  sqrt(  (  xp*onas(l,HP)  -  onas(HP,l)*xp’  )."2  ♦... 

(  yp*onas(l,HP)  -  onao(HP,l)*yp’  ).‘2  ); 
dxpv  ■  -  xp*onas(l,HV)  ♦  onas(HP,i)*xv’ ;  X  dxpv(r,c)  ■  xv_c  -  xp_r  • 
dypv  ■  -  yp*onas(l,HV)  ♦  onas(HP,l)*yv’ ; 
dxvv  •  -  xv*onas(l,HV)  ♦  onas(HV,l)*xv’ ; 
dyw  ■  -  yv*onas(l,HV)  ♦  onas(HV,l)*yv' ; 

App  ■  CFUHpp (Rpp, CLEHGTH ) ;  X  pick  corralation  fun  hara. 

Apu  •  CFUHpu(dxpv, dypv, CLEHGTH) ;  X  pick  corralation  fun  hara. 

Apv  ■  CFUHpv(dxpv, dypv, CLEHGTH) ;  X  pick  corralation  fun  hara. 

Auu  ■  CFUHuu(dxvv, dyw, CLEHGTH)  ;  X  pick  corralation  fun  hara. 

Avv  ■  CFUHvv (dxvv, dyw, CLEHGTH )  ;  X  pick  corralation  fun  hara. 

Auv  ■  CFUHuv (dxvv, dyw, CLEHGTH)  ;  X  pick  corralation  fun  hara. 

A»C  App,  Apu,  Apv; . . . 

Apu', Auu,  Auv;... 

Apv' ,Auv' ,Avv] ; 

L  •  ■P42*MV; 

A(1:L+1:L‘2)  ■  A(1 :L^1 :L‘2) .*(1  ♦  Vn);  Xtaking  cara  of  maasuramant  arrors. 
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dxip  ■  -  zl«on«s(l.lIP)  ♦  oii«a(MI.l)*xp* ;  Xdx(r,c)«xp.c  -  xi.r 

dyip  ■  -  7i*on«s(l,NP)  ♦  on«8(HI,l)*yp’ ; 

dxiv  ■  -  zi*on«8(l,IfV)  ♦  on«a(HI,l)*zv» ;  Xdz(r,c)»zv_c  -  zi.r 

dyiv  ■  -  yi*onea(l,NV)  ♦  onaa(NI,l)«yv' ; 

Cup  >  CFUNup(dzip,dylp,CLENGTH) :  %  pick  corralation  fun  hara. 

Cuu  ■  CFUNuu(dziv,dyiv,CLENGTH) :  %  pick  corralation  fun  hara. 

Cuv  ■  CFUNuT(dziv,dyiT,CLENGTH) ;  %  pick  corralation  fim  hara. 

Cvp  •  CFUNvp(dzip,dyip,CLENGTH) ;  X  pick  corralation  fun  hara. 

Cvu  «  CFUNvu(dziv.dyiv,CLENGTH) ;  X  pick  corralation  fun  hara. 

Cvv  >  CFUNvvCdziv.dyiv.CLENGTH) ;  X  pick  corralation  fun  hara. 

Ci’«[Cup,  Cuu,  Cuv]; 

Cv»[Cvp,  Cvu,  Cvv]; 

X  napping 

ui  ■  Cu*(A\Cp;u;v]); 
vi  ■  Cv*(A\[p;u;v]) ; 

au  *  1  -  diag((Cu/A)*Cu')/CFUNuu(0,CLENGTH) ;  Xpick  corralation  fun  hara. 
av  ■  1  -  diag((Cv/A)*Cv’)/CFUNvv(0,CLEllGTH) ;  Xpick  corralation  fun  hara. 


A. 6  Cross-correlation  Functions 

■■■>  function  y  ■  CFUNppCx.L) 
y  ■  «xp(-(x/L) .‘2) ; 

■■■>  function  Cpu  ■  CFUNpuCdx.dy .CLENGTH) 
lambda  »  l/CLEMGTH‘2: 
r  ■  sqrtC  dx."2  ♦  dy.‘2  ); 

Cpu  ■  dy  .*  (2*lambda  *  CFUNpp(r,CLENGTH)) ; 

■■■>  function  Cpv  ■  CFUNpv(dx,dy,CLENGTH) 
lambda  -  l/CLENGTH‘2: 
r  ■  sqrtC  dx."2  ♦  dy.“2  ); 

Cpv  «  -  dx  .•  (2*lambda  •  CFUllpp(r,CLENGTH)) ; 

■■■>  function  Cup  ■  CFUNup(dx,dy,CLENGTH) 

Cup  *  -CFUHpu(dx,dy,CLENGTH) ; 

«■■>  function  Cuu  ■  CFUIiuu(dx,dy,CLENGTH) 
lambda  »  1/CLEHGTH“2; 
r  ■  8qrt(  dx.'‘2  ♦  dy.“2  ); 

Cuu  ■  (1  -  2*lambda*dy."2)  .*  (2«lambda  •  CFUNppCr.CLEMCTH)) ; 

■»■>  function  Cuv  ■  CFUNuv(dx,dy .CLEHGTH) 
lambda  >  l/CLENGTH‘2; 
r  ■  sqrtC  dx."2  ♦  dy.‘“2  ); 

Cuv  •  2*lambda*dx  .*  dy  .*  (2*lambda  *  CFUNppCr ,CLENGTK)) ; 

■■■>  function  Cvp  ■  CFUIIvp(dx,dy,CLENGTH) 

Cvp  ■  -CFUHpv(dx,dy,CLENGTH) ; 

■»•>  function  Cvu  ■  CFUNvu(dx,dy .CLENGTH) 

Cvu  -  CFUNuv(dx,dy,CLEHGTH'>; 

■■■>  function  Cvv  ■  CFUNvv(dx,dy,CLEHGTH) 
lambda  -  l/CLENGTH‘2; 
r  ■  sqrtC  dx."2  ♦  dy."2  ); 

Cvv  ■  (1  -  2*lambda*dx . '2)  .*  (2*lambda  *  CFUNppCr.CLENGTH)) ; 
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B  Deriving  cross-correlation  functions 

The  relations  in  Table  4  and  5  are  derived  in  this  appendix.  First  recall  the  definitions  of 
7,  Ox  and  Oy". 

7  =  \/(*i  -  *2)^  +  (yi  -  *2  -  *1  =  7  cos  Ox,  1/2  -  yi  =  7  cos  tty. 

From  these  definitions  follow  these  identities: 

^  =  -  cosox;  =  cosox;  ^  =  -  cosa„;  =  cosa,; 


^^7^  = -^(l-cos^ov);  =  icosQxCoso„;  = -i  (l  -  cos^  a,)  • 

By  these  identities  and  Equation  (32),  and  can  be  derived: 

Ft.  =  £:{t«i.,».Wx„»)}  =  =  cco.^. 

Similarly,  we  have 

dF  ,  dF 

^  COS  Qy  OTIC*  ~  cos  Otr  • 

07  07 

The  remaining  four  relations  in  Table  4  follow  from  the  above  two  equations  and  the  iden¬ 


tities  before: 


Fuu  =  ^{«(a;i,yi)t‘(*2.!/2)}  = -£^{^^^^^^«(ar2,J/2)}  = 

dcosoy^  ,,  2 

=  — ^  =  -co8  a„^-(l-cos 


£{u(ii,yi)i>(i2,y3)}  =  -£{ v(i2, ya)}  = 

dcosox^  {<PF  ldF\  ^ 

- ^ 


^{v(*liyi)w(»2.y2)}  =  E  {  «(  »2 . 1/2  )}  = 


dcosoxg  J  <PF  2  ,ldF 

-^:=-co.  a,^-(l-cos’a.)-^. 


and 
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The  results  in  Table  5  are  obtained  by  substituting  F(7)  in  Table  4  by  a  Gaussian 
function: 

/’(7)  =  V'^exp(-A7^ 

Thus  we  have 

F'(7)  =  -2A7F,  F"(7)  =  -(1  -  2A72)  •  2XF 

and  then 

=  7cosay  •  2AF(7)  =  -Fuv»,  F^  =  -7  cos  a*  •  2AF(7)  = 

F„„  =  [cos^ay(l  -  2A7^)  +  sin* Oy]  •  2AF  =  (1  -  2A7*co8*  Oy)  •  2AF(7); 

F„„  =  cos  Ox  cos  Oy  |-(1  -  2A7*)  +  ij  •  2XF  =  2A7*  cos  o*  cos  Oy  •  2AF(7)  =  F^; 

Ftn,  =  [cos*ar(l  -  2A7*)  +  sin*  Qx]  •  2XF  =  (1  -  2A7*cos*  a,)  •  2AF(7). 
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C  Time  series  plots 

Two  sets  of  time  series  plots  are  collected  in  this  appendix.  The  first  set  plots  (Figure  9) 
the  pressure  time  series  which  are  interpolated  to  the  PIES  sites  from  the  final  pressure 
maps  (dashed  lines)  and  the  absolute  pressure  on  the  PIES  sites  determined  in  section  3.3.2 
(solid  lines).  The  second  set  (Figures  10  and  11)  displays  the  time  series  of  geostrophic 
current  on  current  meter  sites  derived  from  the  pressure  maps  (dashed  lines)  and  the  current 
measurements  on  the  corresponding  sites  (solid  lines). 
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Figure  9:  Meuured  presture  (aoUd)  with  mapped  preMure  (daahed) 
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Figure  11:  Measured  (solid)  and  mapped  currents  (dashed)  (v  component) 
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measurements  by  Inverted  Echo  Sounders  (PIES)  and  bottom  current  meter  measurements  in  the 
Central  Array  between  June  1988  and  August  1990. 

A  key  new  feature  is  that  all  dally  maps  can  be  adjusted  to  a  fixed  reference  level  by 
combining  measurements  of  bottom  pressure  and  near-bottom  currents,  (Pjj,  u,v>.  While  objective 
stream  functions  mapped  from  (u,v)  alone  are  spatially  consistent,  the  reference  levels  vary 
temporally  from  map  to  map.  On  the  other  hand,  each  bottom  pressure  record  Pj,  has  a 
temporally  consistent  reference  but  the  reference  can  vary  from  site  to  site.  The  combination 
allows  us  to  adjust  all  to  a  common,  consistent  reference  level  in  multivariate  objective  maps.* 
The  pressure  maps  and  error  maps  are  displayed  daily  for  26  months  with  overlaid  measured 
current  vectors.  Mapping  procedures  are  documented  along  with  error  analyses  and  comparisons 
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19.  ABSTRACT 


with  measured  fields.  The  results  produce  bottom  pressure  and  current  fields  with  typical 
error  of  only  2  mbar  and  2  cm/sec,  compared  to  typical  signal  standard  deviation  of  6  mbar 
and  9  cm/s. 

The  report  is  one  of  a  series  of  data  reports  that  document  aspects  of  our 
participation  in  the  SYNoptic  Ocean  Prediction  (SYNOP)  experiment. 


